Curvature-induced geometric potential in strain-driven nanostructures 
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We derive the effective dimensionally reduced Schodinger equation for electrons in strain-driven 
curved nanostructures by adiabatic separation of fast and slow quantum degrees of freedom. The 
emergent strain-induced geometric potential strongly renormalizes the purely quantum curvature- 
induced potential and enhances the effects of curvature by several orders of magnitude. Applying 
this analysis to nanocorrugated films shows that this curvature-induced potential leads to strongly 
enhanced electron localization and the opening of substantial band gaps. 

PACS numbers: 02.40.-k, 03.65.-w, 68.65.-k, 73.22.-f 



The experimental progress in synthesizing low- 
dimensional nanostructures with curved geometries - the 
next generation nanodevices [IH3] - has triggered the in- 
terest in the theory of quantum physics on bent two- 
dimensional manifolds. The current theoretical paradigm 
relies on a thin-wall quantization method originally intro- 
duced by Da Costa 4 . It treats the quantum motion on a 
curved two-dimensional (2D) surface as the limiting case 
of a particle in three-dimensional (3D) space subject to 
lateral quantum confinement. With Da Costa's method 
the surface curvature is eliminated from the Schrodinger 
equation at the expense of adding a potential term to 
it. This simplifies the problem substantially, as elec- 
trons now effectively live in 2D space, in presence of a 
curvature-induced potential. 

As the curvature induced potential is entirely of quan- 
tum origin - it is in magnitude proportional to h - its 
physical consequences in condensed matter systems can 
only be observed on the nano-scale. In this realm the 
curvature-induced geometric potential can cause intrigu- 
ing phenomena, such as winding-generated bound states 
in rolled- up nanotubes [5j |6] and topological bandgaps in 
periodic curved surfaces [7]. However, the experimental 
realization and exploitation of such phenomena is hin- 
dered by the fact that in actual systems with curvature 
radii of hundreds of nanometers, Da Costa's quantum 
potential is still very weak and typically only comes into 
play on the sub-kelvin energy scale. 

In this Letter, we mitigate this problem by developing 
a thin-wall quantization procedure which explicitly ac- 
counts for the effect of the deformation potentials of the 
model-solid theory . By employing a method of adia- 
batic separation of fast and slow quantum degrees of free- 
dom, we show that local variations of the strain render an 
emergent geometric potential that strongly (often gigan- 
tically) boosts the purely quantum geometrical potential. 
The theoretical framework that we develop is immedi- 
ately relevant for electronic nanodevices as the present- 
day nanostructuring method [T] is based on the tendency 
displayed by thin films detached from their substrates to 
assume a shape yielding the lowest possible elastic energy. 
As a result, thin films can either roll up into tubes pUlfTT] 



or undergo wrinkling to form nanocorrugated structures 
[21 [121 [13] . A key property of such bent nanostructures is 
the nanoscale variation of the strain. It for instance leads 
to considerable band-edge shifts [1_4^ with regions under 
tensile and compressive strain shifting in opposite direc- 
tion. To include such a nanoscale sequence of potential 
wells in the elastically relaxed structure one obviously 
needs to go beyond Da Costa's thin-wall quantization 
framework since it intrinsically couples the transversal 
quantum degrees of freedom to the tangential quantum 
motion along the curved surface. 

In the thin-wall approach quantum excitation energies 
in the normal direction are raised far beyond those in the 
tangential direction by the lateral confinement. Hence 
one can safely neglect the quantum motion in the nor- 
mal direction thereby reaching an effective dimensionally 
reduced Schrodinger equation. As opposed to a classical 
particle, a quantum particle constrained to a curved sur- 
face retains some knowledge of the surrounding 3D space. 
In spite of the absence of interactions, it indeed experi- 
ences an attractive potential of geometrical nature [1]. 
It has been shown that Da Costa's thin-wall quantiza- 
tion procedure is well-founded, also in presence of exter- 
nally applied electric and magnetic fields |T5] . Even more 
the experimental realization of an optical analogue of the 
curvature-induced geometric potential [16| has provided 
empirical evidence for the validity of this approach. 

To develop a theoretical framework for strain-induced 
geometric potentials we therefore use the same concep- 
tual framework. We start the mathematical description 
by defining a 3D curvilinear coordinate system [see Fig. [l] 
for a generic bent nanostructure. The stress-free surface 
S is parametrized as r = r(s,y) where y is the coor- 
dinate along the translational invariant direction of the 
thin film and s is the arclength along the curved direction 
of the surface S measured from an arbitrary reference 
point. Nothing prevents the thin film to be deflected 
also along the y direction. However, in the remainder 
we will neglect this effect for simplicity. The 3D por- 
tion of space of the thin film can be then parametrized 
as R(s, y, q^) — r(s, y) + q3N{s) with N{s) the unit vec- 
tor normal to S. We can evaluate the strain distribution 
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FIG. 1: Schematics of the three-dimensional coordinate sys- 
tem for a strain-driven bent nanostructure with positive ra- 
dius of curvature R and total thickness S. The qs — surface 
5 corresponds to the stress-free surface. Regions with qs > 
{qs < 0) are under compressive (tensile) strain. 



in the thin film by assuming that along the undeflected 
direction Cy = 0. It is well known Tf! that the strain 
in the direction along the surface varies linearly across 
the thin film as Cs = —q^K{s) where k{s) is the princi- 
pal curvature of the surface S. The strain in the normal 
direction can be related to eg by means of the Poisson 
relation eq^ = — [v/ (1 — v)] with v the Poisson ra- 
tio. From the linear deformation potential theory |8] we 
then have that the strain induced shift of the conduction 
band corresponds to a local potential for the conducting 
electrons Ve( 5,(73) = 7'3'3k(s) with 7 > yielding an at- 
traction toward regions under tensile strain [c.f. Fig. [l]. 
The characteristic energy scale 7 which can be explic- 
itly computed in the different conduction valleys of the 
nanostructure is proportional to the shear and the hydro- 
static deformation potentials and typically lies in the eV 
scale for conventional semiconductors By adopting 
Einstein summation convention, the Schrodinger equa- 
tion for the quantum carriers in the effective mass ap- 
proximation then takes the following compact form |15) 
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In the equation above G'-' corresponds to the 3D metric 
tensor of our coordinate system, the covariant derivative 
Di is defined as "Di = diVj — ^ijVk with vj the covariant 
components of a generic 3D vector field, and the P^j 's 
are the Christoffel symbols. By expanding Eq. ([T]) by 
covariant calculus 15 we get 
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where we defined Hs = 1 — k(s) q'3 , and we introduced 
a squeezing potential in the normal direction Vxi^q^). In 
the following it will be considered as given by two infi- 
nite step potential barriers at ±(5/2 where 5 is the total 
thickness of the thin film. However, our results can be 



straightforwardly generalized to other types of squeezing 
potential, e.g. harmonic traps. 

In the same spirit of the thin-wall quantization pro- 
cedure [1] we next introduce a new wavefunction x 
for which the surface density probability is defined as 
/ Ixls, 93)prf93- Conservation of the norm requires 

1 /2 

X = -0 X Hg . The resulting Schrodinger equation is 
then determined by the Hamiltonian 
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Assuming the thickness of the thin film 5 to be small com- 
pared to the local radius of curvature R{s) — k(s)~^, the 



Hamiltonian Eq. (|3| can be expanded as "H = J2k Is 
At the zeroth order in q^ we recover precisely the effec- 
tive Hamiltonian originally introduced by Da Costa [4] 
which disregards the strain-induced shifts of the conduc- 
tion band and guarantees the separability of the tangen- 
tial motion from the transverse one. In this case the effect 
of the curvature eventually results in the well-known geo- 
metric potential Vg(s) = — /i^k(s)^/(8 to*). Strain effects 
can be explicitly monitored by retaining linear terms in 
93 in which case we obtain the effective Hamiltonian 
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The strong size quantization along the normal direction 
allows us to employ the adiabatic approximation and 
solve the Schrodinger equation for the effective Hamil- 
tonian Eq. ^ considering the ansatz for the wavefunc- 
tion xl^,?/, (73) — X^i^jQs) ^ X^i^iU) where the normal 
wavefunction x^ solves at fixed s the one-dimensional 
Schrodinger equation for the "fast" normal quantum de- 
grees of freedom 
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^ (5) 
Here i indicates the transversal subband index and 
7(5) = -f — h"^ [i^isY + 9sK{s)/k{s)] /(4to*) corresponds 
to the typical energy scale of the deformation potential lo- 
cally renormalized by curvature effects. The latter term 
yields an attraction towards regions under compressive 
strain (where the local curvature is higher) competing 
with the dominant strain-induced attraction towards the 
tensile regions of the nanostructure [c.f. Fig. |2ja)]. 

The Hamiltonian for the slow tangential quantum mo- 
tion can be found by first integrating out the 93 quan- 
tum degree of freedom and then performing the addi- 
tional rescaling of the tangential wavefunction ~^ 
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FIG. 2: (a) Schematics of the electron density probabiUty 
in the lowest normal eigenstate as a function of the coordi- 
nate gs measured in units of the thin film thickness S for a 
positive value of the local curvature in a strain-free (7 = 0) 
and a strained (7 7^ 0) curved nanostructure. The vertical 
dashed line indicates the stress-free surface 53 = 0. (b) De- 
viation of the locally renormalized effective mass Am^(s) = 
(m|;(s) — m*)/m* as a function of the typical energy scale 
7(3) for different values of the expansion parameter k(s)5 . 
(c) Behavior of the adiabatic potentials (s) for the first 
two transversal subband measured from the quantum well 
levels ft^7r^i^/(2m*(5^) for k{s)5 = lO"'^ as a function of 7(3). 
The continuous lines indicate the exact behavior whereas the 
points are the results of the perturbation theory approxima- 
tion. 



X ^yi + 2K{s){x^\q3\x^). The final form of the di- 
mensionally reduced tangential Hamiltonian is as follows 
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where the . . . indicate the diagonal adiabatic corrections 
and we introduced the locally renormalized effective mass 



mjj(s) = m7[l-H2«:(s)(x^|g3|x'^)]- Fig. gb) shows 
the deviation of the locally renormalized effective mass 
from its bare value for different values of the expansion 
parameter k{s)5. The strain-induced localization of the 
normal wavefunction in the tensile regions of the nanos- 
tructure leads to an heavier effective local mass enhanced 
at the points of maximum curvature. This is at odds 
with the lighter local effective mass one would find in 
the absence of strain effects. Even for values of the lo- 
cal energy scale 7(5) much larger than the characteristic 
energy of the normal quantum well energy h?' / {2m* 5"^) , 
this local renormalization of the effective mass is so small 
that for all practical purposes the use of the bare effec- 
tive mass in Eq. (|6| is justified. More rigorously one can 
show the reliability of this approximation in the regime 
7(s)k(s)5 < It? / {2m* 5"^) where the strain-induced lin- 
ear potential appearing in the fast Schrodinger equation 
Eq. ^ can be treated perturbatively. Since such a condi- 
tion is typically satisfied in conventional semiconducting 
nanostructures with a total thickness in the nanometer 
scale, we will limit ourselves to this regime from here 
onwards. 

In Fig. [2]^c) we show the behavior of the adiabatic 
potentials in the first two transversal subband mea- 
sured from the normal quantum well levels i??^'s. The 
characteristic quadratic dependence on the local energy 
scale 7(5) is well reproduced by the analytical formula 
Ef{s) = E^^ + 2m* S'^j{s)^/h^ x fi ^(s)^ with the nu- 
merical constants that can be calculated using second- 
order perturbation theory as /i ~ — 10~^, /2 ~ 3 • 10~*, 
etc. From this it also follows that the distance among the 
potential energy surfaces E^ {s) - E^ {s) > E^^ ~ E^^ 
thereby guaranteeing the reliability of the adiabatic ap- 
proximation in the regime k{s)5 ^1. As a result, we 
then find an effective Hamiltonian for the electronic mo- 
tion along the curved nanostructure which in the first 
transversal subband reads 
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where we left out the constant energy term /i^tt^/ {2m*6^) 
and we neglected the diagonal adiabatic corrections. The 
second term in Eq. [7] corresponds to a geometric poten- 
tial whose strength is renormalized by strain effects as 
VR = 1 + 4|/i| {2m*5'^-i/h^f . Remarkably we find this 
renormalization to become extremely large in case of the 
natural hierarchy of energy scales 
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Considering for instance a value of 7 ^ 10 eV, a charac- 
teristic quantum well energy h?/{2m*6^) ^ 5 meV and 
the typical tangential kinetic energy fi^ / [2m* R{s)'^] ^ 
1/ieV, we find an enhancement of the curvature-induced 
potential by ^ 10*. As we show below, this gigan- 
tic renormalization of the geometric potential has pro- 
found consequences on the electronic properties of low- 
dimensional nanostructures with curved geometry. It is 
worth noting that in absence of strain effects (7 — >■ 0), 
Eq. ([7]) corresponds to the effective tangential Hamilto- 
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FIG. 3: (a) Sketch of the stress-free surface of a nanocorru- 
gated thin film, (b) Bandstructure of a nanocorrugated film 
with Aq = 0.2 and different values of the geometric potential 
renormalization factor vr. We show only the first two bands 
close to the bottom of the conduction band, (c) Behavior of 
the topological bandgap as a function of the dimensionless pa- 
rameter Aq for different values of vr. (d) Tangential ground 
state density probability for Aq = 0.2 and different values of 
the renormalization strength of the geometric potential. 



nian introduced by Da Costa augmented with an higher- 
order curvature induced geometric potential cx k(s)^ aris- 
ing as a consequence of the finite thickness of the thin- 
film nanostructure. 

We now use this theoretical framework to analyze the 
influence of the strain-induced geometric potential on the 
electronic states of a nanocorrugated thin film with pe- 
riod 27r/q and total thickness 5 [c.f. Fig. [sjja)]. The 
stress-free surface can be parametrized in the Monge 
gauge as 53 {x) = A cos q x where A is the amplitude 
of the corrugation. In the shallow deformation limit 
Aq <C 1, we can express the arclength of the layer s ~ a; 
whereas the local curvature of the stress free surface 
k(s) ~ —A <f cos s). Thus, the problem reduces to the 
"flat" motion of free electrons embedded in a curvature 
induced periodic potential. It is straight-forward to ob- 
tain the energy spectrum of the "slow" tangential motion 
of Hamiltonian Eq. ([t]) in the first transversal subband 
for a thin film thickness much smaller than the corru- 
gation wavelength as in this case the last two terms in 
Eq. ^ can be neglected. Fig. [sj^b) shows the behavior 
of the first two bands for different values of vji. The zero 
of the energy has been chosen as the bottom the conduc- 
tion band. The qualitative behavior of the bandstructure 
is maintained when strain effects are taken into account, 
but now with curvature induced gaps |19j at momenta 
k — ±q. The gaps increase quadratically in magnitude 
with the energy scale of the deformation potential 7. In- 
deed we find AE ~ vrH'^A'^ q'^/^lQm*), in agreement 
with the numerical analysis [c.f. Fig. [sjjc)]. In Fig. [sjjd) 
we show the ground state density probability for differ- 



ent values of vr. By increasing strain effects, one finds 
a continuous crossover from extended-like states to elec- 
tronic states localized precisely at the points of maximum 
and minimum curvature, which is in accordance with the 
results of a purely numerical approach [T5] . 

In conclusion, by employing a method of adiabatic sep- 
aration of fast and slow quantum degrees of freedom, we 
have derived a dimensionally reduced Schrodinger equa- 
tion in strain-driven nanostructures. The strain effects 
render an often gigantic renormalization of the curvature 
induced geometric potential which has profound conse- 
quences upon the electronic properties of these materi- 
als. Applying our theoretical framework to the case of 
nanocorrugated thin films we find an enhanced electron 
localization and the opening of substantial band gaps on 
an experimentally relevant energy scale. It can also be 
applied to, for instance, 2D nanotubes roUed-up in the 
shape of an Archimedean spiral, where the effect of the 
geometric quantum potential leads to bound states whose 
number coincides with the winding number [5]. The in- 
clusion of strain effects will lead to a proliferation of such 
bound states, strongly affecting the electronic and trans- 
port properties of these roUed-up nanostructures. 
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